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Abstract 

We investigate the non-perturbative results of multi-dimensional forced Burg- 
ers equation coupled to the continuity equation. In the inviscid limit, we derive 
the exact exponents of two-point density correlation functions in the universal 
region in arbitrary dimensions. We then find the universal generating function 
and the tails of the probability density function (PDF) for the longitudinal veloc- 
ity difference. Our results exhibit that in the inviscid limit, density fluctuations 
affect the master equation of the generating function in such a way that we can 
get a positive PDF with the well-known exponential tail. The exponent of the 
algebraic tail is derived to be —5/2 in any dimension. Finally we observe that 
various forcing spectrums do not alter the power law behaviour of the algebraic 
tail in these dimensions, due to a relation between forcing correlator exponent 
and the exponent of the two-point density correlation function. 
PACS numbers 47.27.Gs and 47.40.Ki 

The interest in solving the randomly driven Burgers equation with a large scale 
driving force is motivated by the hope that it can provide us with the first solv- 
able model of turbulence. Consequently tremendous activity emerged on the non- 
perturbative understanding of Burgers turbulence |2|, |3|, ^ ||, ^, ^, ^, 10, 11, 12 1. 



As an original attempt, Polyakov offered a field theoretical method to derive the prob- 
ability distribution or density function (PDF) of velocity difference in the problem 
of randomly driven Burgers equation in one dimension [|l|. The problem of compu- 
tation of correlation functions in the inertial range is reduced to the solution of a 
closed partial differential equation ||l|, However, Polyakov's approach was based 
on the conjecture of the existence of the operator product expansion (OPE). This 
method was then extended to the forced Burgers equation coupled to the continuity 
equation and some results were derived for the behaviour of the probability density 
function tails and for the value of intermittency and density-density correlators expo- 
nents [S, Ol. On the other hand, some extensive numerical simulations show that the 
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predictions of these theoretical works coincide with the numerical simulations, within 
a good approximation |10, I l3|. 

|, |§|, 0, |, H, [l^, m that the PDF in Burgers turbulence 



It was known ||ll ^, 



behaves physically different in two different intervals of velocity difference u. When 
\u\> Urms the PDF behaves as: 

P(u,r)=rG(-^) 

'-^rms 

which depends on the single-point property Urms (root-mean-square of velocity fluc- 
tuations) and therefore is not a universal function. While in the interval u -C Urms 
and r <^ L (where L is the energy input scale), the PDF can be represented in the 
universal scaling form: 

^/ X 1 ^, n , 

where F{x)dx = 1 and the exponent 5 is related to the random force spectrum 
[p!^]. Results for one-dimensional case indicates that the function F{x) behaves like 
exp[— c((5)x^] when x — > -|-oo, where the coefficient c{5) can be evaluated from the 
theory |^ In this region, various algebraic behaviours were predicted for the 
PDF in the limit x — > — oo and determination of the asymptotic behaviour of the 
PDF is at the moment controversial. Several different proposals have been made, 
each leads to an asymptotic experssion of the forms ~ \u\~^ but with a variety of 
values for the exponent a includes 2 |, 0], the range 5/2 to 3 0, |, |l2|, 3 |0| 
and 7/2 |@, |8|, |9|. Numerical simulations performed in ||l^ shows that the power law 
exponent of the PDF depends on the forcing spectrum. The investigations of forced 
Burgers turbulence models have also given furthur understanding of intermittency in 
the velocity structure functions. It is believed that intermittency in these systems is 
a consequence of the algebraic tail of PDF implying the scaling exponent of velocity 
structure functions, is saturated to the value ^„ = 1 for n > nc, where the value 
of critical moment number Uc depends on the forcing function spectrum This 
simply means that the moments Snyuc dominated by the non-universal region 
and thus vary with the sigle-point property Urms induced by the large scale random 
force. At the same time, the moments with n < ric are universal [ffl, O]. 



The first attempt on d-dimensional Navier-Stokes turbulence was made in |15] as 
a renormalized perturbation expansion with l/d as a small parameter. This attempt 
failed because 1/d factors appearing in the expansion are cancelled by the 0{d) mul- 
tipliers resulting from the summation over the d components of velocity field. The 



problem of d-dimensional turbulence was recently revisited |16] using Polyakov's non- 
perturbative approach. It was shown that in the limit d — > oo Kolmogorov scaling 
is an exact solution of the incompressible Navier-Stokes equations and the pressure 
contributions to the equations are responsible for a clear distinction between the 
Navier-Stokes and Burgers dynamics. In an other work the asymptotics of the 
PDF of velocity gradient were found in d-dimensional randomly driven Burgers equa- 
tion in a compressible fluid. 
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In the present paper we investigate the isotropic forced Burgers equation coupled 
to the continuity equation in arbitrary dimensions. In the inviscid hmit, we find a 
complete closed equation for the generating function of velocity difference. We solve 
this equation for the longitudinal component of velocity difference in the universal 
region and show that the main requirments on the PDF fix the exponent of two- 
point density correlation function. This results in a dimension dependent exponent 
for density fiuctuations. This then means that the equation governing the PDF is 
independent of dimensions. In this region, our results predict the exponential decay- 
ing of the right tail of PDF of velocity difference which is also obtained in several 
other works |, |, |, |, |, g |, 0, ||]. Thus the power law decay of the left tail is 
independent of the dimensionality and the forcing spectrum with the exponent —5/2 
in any dimension. While the previous results for one-dimensional forced Burgers tur- 
bulence exhibit a force spectrum dependent tail [Q, ^, 13| or a different power law 
decaying 0, ^, ^ Finally, we discuss Polyakov's approach for dealing with 
the problem in the limit of infinitesimal viscosity. For a special forcing spectrum 
this method would not give anything more than the inviscid calculations. Our work 
extends a previous result |12] on two and three-dimensional Burgers turbulence to 
higher dimensions. 



Let us start with the Burgers and continuity equations: 

9tu + (u-V)u = i/V2u + f(x,t) (1) 

dtp + V ■ (pu) = (2) 

for the irrotational velocity field u(x, t), viscosity u and density p m. d dimensions. 
The force f (x, t) is the external stirring force, which injects energy into the system on 
a length scale L. More specifically, one can take, for instance a Gaussian distributed 
random force, which is identified by its two moments: 



(/m(x>*))=0 

(/^(x, t)/,(x', t')) = 6{t - t')k^,ix - x') (3) 

where p,u = 1,2, ... ,d and the correlation function k^^(r) is concentrated at some 
large scale L. The problem is to understand the statistical properties of the velocity 
and density fields which are the solutions of equations (1) and (2). Following Polyakov, 
we consider the following two-point generating functional: 



F2(Ai, A2,xi,X2,t) = (/)(xi,t)/9(x2,t)exp(Ai • u(xi,t) + A2 • u(x2,t))) (4) 

where the symbol (. . .) means an average over various realization of the random force. 
Now one can show that F2 satisfies the following equation: 
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dtF2 + 



d 



i=l,2;M=l,...,rf ^■^'^'^ 



dn,F2 - 



E 

*J=l,2;/i,i^=l, 



(5) 



where the first two terms on the left-hand side of equation (5) come from the terms on 
the left-hand side of equations (1) and (2), and the third is the contribution of forcing 
term, in which we have used Furutsu-Novikov-Donsker formula |17]. Also, L'2-term is 
the contribution of dissipation. D2 is the anomaly term and has the following form: 



D2 = {up{xi)p{x2)[Xi ■ V2u(xi) + A2 • V2u(x2)] exp(Ai • u(xi) + A2 • u(x2))) (6) 

It is noted that although the advection contribution are accurately accounted for 
in this equation, it is not closed due to the dissipation term. In what follows, we 
consider the inviscid (v = 0) Burgers equation to avoid the anomaly problem and 
leave Polyakov's approach for dealing with the problem in the limit 1/ — > to the last 
part of the present paper. Now we change the variables as xj- = xi ibx2, A+ = Ai -|- A2 
and A_ = By the Galilean invariance and the spatial homogeneity assumptions 

the variables x+ and A+ can be set to zero. It is found from the equation (5) that in 
the stationary state we have: 

^ 9^5^^' + ^ ^ A,.A.[V(x)-V(0)]i^2 = (7) 

li=l,...,d ^ ^^ ii,i/=l,...,d 

where we have used x and A instead of x_ and A_ for simplification. Also, because of 
isotropy F2 can depend only on the absolute values of vectors x and A and the angle 9 
between them as F2 = F2{r, A, s) where r = |x|, A = |A| and s = cosO = Q2 Xfj_\p)/r\. 
We suppose the stirring correlation function has the following form: 

= - ^^'(V + (8) 

with fci and L = 1. Now, using spherical coordinates (r. A, s) it can be shown from 
equations (7) and (8) that F2 satisfies the following closed equation for homogeneous 
and isotropic turbulence: 



[sdrdx r d; H r ds H : drds 

rX rX X 

+ ^^~' h xds-r^X^{l + 2s^)]F2 = Q (9) 
r 

The one-dimensional case of equation (7) is easily recovered by setting s = 1. We wish 
to consider the longitudinal velocity component statistics which results from equation 
(9) by taking the limit s — > 1. Assuming F2 be a regular function near s = 1 we can 
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safely drop the terms multiply by (1 — s^) in equation (9). We propose the universal 
scale invariant solution of equation (9) in the following form: 



F2ir,X,s)=g{r)F{Xr\s) 

g{r) = r""-* (10) 

In equation (10), g{r) is the conditional two-point correlation function of density 
field conditioned on a fixed value of velocity differnce. We assume that g{r) depends 
only on r and its dependence on the velocity interval appears in the exponent a^, 
we shall discuss this further in the next section. For a general stirring correlation 
function k^u^{^ — f^): the exponent 5 is found by substituting the proposed form of 
generating function as 5 = which in our case (C = 2) is = 1. Indeed, we assume 
that the two-point density correlation function exists, and therefore it is necessary 
to find F{Xr^,s) such that it tends to a constant in the limit of A — > 0. Also, It is 
straight forward to show that the mean value of velocity difference is zero as expected 
from the homogeneity and isotropy constraints. Now let us consider the longitudinal 
velocity component statistics suggested by equation (9) in the limit s — > 1. In this 
limit, we assume that the generating function of velocity difference F{Xr, s) has the 
following form: 



F{Xr,s) = F{Xrs) (11) 

this form ensures the factorizing property of the angular part of structure functions 
as Sn{r,s) ~ s"S'„(r) when n < 1. The factorization of the angular part of velocity 
structure functions in the limit s ^ 1 has also been known for the Navier-Stokes 
turbulence jl^. Plugging the ansatz (10) and (11) into equation (9) and rewriting it 
for the variable z = Xrs gives the following equation for the generating function of 
longitudinal velocity difference: 



zF"{z) + {d- ad)F'{z) - 'iz^F{z) = (12) 

It is interesting that equation (12) is similar to the equation first derived by Polyakov 
for the problem of one dimensional Burgers. In that work the effect of the viscous 
term is found by taking the limit of infinitesimal viscosity applying the self-consistent 
conjectures of operator product expansion. The anomaly terms which arise from the 
viscous term modify the master equation of the generating function in such a way 
that a positive, finite and normalizable PDF can be found. Comparing equation 
(12) with Polyakov's result shows that the anomaly term in Polyakov's approach 
is replaced with the exponent of two-point density correlation function such that a 
simple change of the parameters maps equation (12) to one derived by Polyakov Q. 



Therefore, as mentioned in |12], the presented approach for inviscid forced Burgers 
turbulence shows that considering the density fluctuations coupled to the velocity 
field, alters the governing equations in such a way that we can obtain positive PDF 
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even in the inviscid hmit and with no need for the viscosity anomaly. It is easy to 
show using equation (12), as discussed in [|l|, that the main requirments on the PDF 
forces us to take the two point density correlator exponents in arbitrary dimensions 
as: 



, 1 / ^ 

ad = d+- (13) 

which yields F(z) = exp[;^(z^/^)]. This result consequently gives the right tail of the 

PDF of longitudinal velocity differrence as ^ exp[— c(^)^] (for ^ +oo) where c = ^. 
This tail coincides with the result of several other approaches [0, ^, ^, ^ P, R, ^, 10|. 
The left tail of the PDF is obtained in the hmit of ^ ^ — oo as or, by 

equation (13), |n|~^/^. As a result, density fluctuations which couple to the velocity 
field appear themselves in the exponent of left tail of PDF and lead to a dimensionality 
independent decay law. 



Now we wish to focus our attention on the dependnce of these results on the 
stirring force spectrum. We consider a general large-scale stirring force correlator as 
follows: 

V(x) = ko6^u - ^rf(V + 2^) (14) 

With this choice the velocity field remains irrotational also the coefficients in expres- 
sion (14) are taken such it simply reduces to equation (8) when ^ = 2 . The last 
term in the equation (5) can be written in the form —ar''X'^F2 for a general form of 
the stirring force. The coefficient a depends on the parameters in the expression of 
the stirring correlator. For the choice as in equation (14) a is (1 -|- 2s^). However the 
exact form of a does not affect the main features of the problem such as algebraic tail 
of the PDF. Taking the universal solution (10) for the new master equation gives the 
following eqation for the generating function F{Xr^s) = F{z) in the limit s — > 1: 



6zF"{z) + {d-ad + S- l)F'{z) - 3z^F{z) = (15) 

where 5 = The PDF of longitudinal velocity difference is the inverse Laplace 

transform of F{z). The positivity, finitness and normalizability requirements on the 
PDF will fix the exponent of two-point density correlator as: 



arf = d+i(C-l) (16) 

also we obtain F{z) = exp[-|j(z3/2)]. The right tail of the PDF can be readily 

deduced as exp[— cx^] when x = ^ -|-oo and c = — |. On the other hand, the limit 
j; — > — oo gives the left tail in the following form: 
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P(n,r) ~ |nr("''-'^+^+i)/^ (17) 

Substituting the value of ad from equation (16) and 5 = we obtain the algebraic 
tail as: 



P{u,r) ^r^^/'^\u\-^/^ (18) 

We observe that the density fluctuations affect the algebraic tail of the PDF and 
modify its decaying exponent such that it becomes independent of the dimension and 
the forcing spectrum. 



As mentioned previously, in the limit — > 0, Polyakov formulated a method for 
analyzing the inertial range statistics based on the conjecture of the existence of OPE. 
Extending the assumptions of OPE to take into account the anomaly term in the case 
of arbitrary dimensions for Burgers equation coupled to the continuity equation, it 
was shown that D2-term has the following structure [Q: 

D2 = aF2 (19) 

where a is generally a function of Ai and A2. Therefore keeping the viscosity infinites- 
imal but nonzero produces a finite effect and a new term on the right hand side of 
equation (9) as: 

[sdrdx r H r ds H r drds 

rX r\ A 

+ ^-^—^dxds - r^\\l + 2s^)]F2 = a{X)F2 (20) 
r 

The A dependence of a (A) anomaly must be chosen to conform the scaling and can be 
changed depending on the properties of the force correlation function. For a general 
correlation function /c^i,~(l — r'') the A dependence of a is fixed as ooA"" where a = 
and oq is a constant. It is evident that in the case of C = 2, a(A) is independent of 
A. In this case, one can easily show from the master equation that the parameter ao 
depends linearly on the mean value of velocity difference field of flow and therefore 
vanishes in the homogeneous isotropic turbulence. However for different types of 
correlations of the stirring force, e.g.A:^j^~(l — r^) with ^ 2, we have to assume 
non-zero ao Q. Thus, for the stirring correlation as the type 1 — r^, the results of 
Polyakov's formalism (in the limit of infinitesimal viscosity) will not be different from 
the inviscid results of homogeneous isotropic turbulence. 

To our knowledge, there is a little informations about the statistics of the density 



field. There exists one simulation [18| to find the PDF tails of density field for one 
dimensional decaying Burgers turbulence in the zero viscosity limit, which exhibits 
a power law tail for the density PDF in the high density regime. In one dimension. 
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Boldyrev |^] has reported a simulation in which the exponent of two-point density 
correlation function is predicted to be ^ 2. 



We would like to thank J. Davoudi, S. Moghimi-Araghi and M.R. Rahimi Tabar 
for useful discussions. 
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